Exact Scale-Invariant Background of Gravitational Waves from Cosmic Defects 
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We demonstrate that any scahng source in the radiation era produces a background of gravita- 
tional waves with an exact scale-invariant power spectrum. Cosmic defects, created after a phase 
transition in the early Universe, are such a scaling source. We emphasise that the result is inde- 
pendent of the topology of the cosmic defects, the order of phase transition, and the nature of the 
symmetry broken, global or gauged. As an example, using large-scale numerical simulations, we 
calculate the scale invariant gravitational wave power spectrum generated by the dynamics of a 
global O(A^) scalar theory. The result approaches the large A'' theoretical prediction as A*""^, albeit 
with a large coefficient. The signal from global cosmic strings is 0(100) times larger than the large 
A^ prediction. 

PACS numbers: 98.80.-k, 98.80.Cq, 04.30.-w, 11.27.-|-d 
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Information about new physics in the very early uni- 
verse is directly accessible through weakly-interacting 
relics. The weakest of all interactions is gravity, and suf- 
ficiently energetic processes leave behind characteristic 
signatures in relic gravitational waves (GW). New GW 
observatories are under construction or have been pro- 
posed These are very likely to detect astrophysical 
sources soon, and have the potential to detect cosmologi- 
cal sources or strongly constrain early universe scenarios. 

Gravitational waves are produced whenever there is 
an energy-momentum tensor with a transverse-traceless 
(TT) part [J]. In the absence of any source, GW are 
also generated quantum mechanically during inflation Q . 
In this paper we study the GW background generated 
by a scaling source, whose energy-momentum tensor is 
proportional to the square of the Hubble parameter. 

Scaling is exhibited by cosmic defects, products of a 
phase transition in the early universe If the vacuum 
manifold M is topologically non-trivial, i.e. has a non- 
trivial homotopy group 7r„(A1) 7^ I, topological field con- 
figurations arise: strings for n= 1, monopoles for 71 = 2, 
and textures for n = 3. For higher n, non-topological field 
configurations are created. When the spontaneously bro- 
ken symmetry is global, the defects generated are termed 
global. When the broken symmetry is gauged, local 
defects appear. Global defects (independently of their 
topology), as well as cosmic strings (global or gauged) ex- 
hibit scaling behaviour, sufficiently long after the phase 
transition that created them |5l-[l(|. All cases, topologi- 
cal or not, local or global, are termed cosmic defects. 

Based on dimensional grounds and causality, it was ar- 
gued in JTj that global phase transitions generate an ap- 
proximately scale-invariant GW background. The ampli- 
tude was estimated with the quadrupole approximation 
for GWs, without any reference to the number N of com- 
ponents of the symmetry breaking field. In the context 
of the large N limit of a global phase transition 0, Q and 



using a full treatment of the tensor metric perturbation 
representing GWs (i.e. not resorting to the quadrupole 
approximation), it was demonstrated in (T^. [isj that in- 
deed an exact scale-invariant background of GW is gen- 
erated by the self-ordering process of the non-topological 
global textures arising after the phase transition. 

In this letter we generalize these results. We demon- 
strate that any scaling source in the radiation era pro- 
duces a background of GWs with a scale- invariant energy 
density power spectrum. In the case of defects, we em- 
phasise that the result is not related to their topology, 
the order of phase transition or the global/local nature of 
the symmetry-breaking process; it is just a consequence 
of scaling. As an example, using lattice simulations as an 
input, we calculate numerically the GW amplitude from 
a system of global 0{N) defects. We consider simula- 
tions with A^ = 2, 3, 4, 8, 12 and 20, thus probing both 
topological and non-topological global defect scenarios. 
We provide evidence that the numerical results converge 
to the large A^ calculation [Eq. ([T|) below] as N~'^, al- 
beit with a large coefficient. The GW power from global 
strings (A^ = 2) is significantly above the trend, around 
100 times greater than the large A^ prediction. 

On the numerical side, our present work is comple- 
mentary to that recently presented in p^ . where nu- 
merical simulations are also performed and qualitative 
agreement with the large A^ result is obtained. There, 
the authors conclude that global defects created after a 
second-order phase transition generate a scale-invariant 
GW background. Such a conclusion, however correct, 
is not fully supported by their numerical spectra: their 
GW spectra could be consistent with scale invariance, 
but there are large fiuctuations. Our aim in this letter is 
precisely to demonstrate that scale-invariance is indeed 
exact for any scaling source, and that this is not related 
to the type of phase transition. 

In the global O(A^) theory, although the field equations 
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are non-linear, analytic calculations are possible in the 
limit of large iV Q , which show that the self-ordering 
dynamics of the non-topological defects exhibits scaling. 
It is also possible to calculate analytically the GW power 
spectrum [ll|, (see also [271]), as 



(1) 



where pc is the critical energy density today, and Mp « 
1.22 X 10^® GeV is the Planck mass. The parameter v 
is the vacuum expectation value (vev) of the scalar field 
and r^rad — 4 X 10~^ is the radiation-to-critical energy 
density ratio today. The amplitude of this background 
does not depend on the frequency /, it is an exact scale- 
invariant background. There is no dependence either on 
the self-coupling A of the symmetry-breaking field. This 
is because the effective theory of the Goldstone modes, 
responsible for the creation of the GWs, is a non-linear 
cr-model, and the coupling disappears when the scalar 
field mode is integrated out. 

Important questions are raised, which we address in 
this letter. How does the scale-invariant GW spectrum 
come about? How does the GW spectrum look in the 
case of topological defects? How does the true GW signal 
from global non-topological textures approach the large 
N result? The last two questions are particularly relevant 
in string-inspired models such as ^15. . Jj6,] , which can have 
(approximate) global symmetries with low N. 

In the following we shall assume that the total energy- 
momentum tensor T^i^ has contributions from ideal ra- 
diation, matter, and defects, so that we can split it as 
T^^ = T^ad _^ ymat _^ ydef ^ Unlike radiation and matter, 
the energy-momentum tensor of defects is not a perfect 
fluid, and supports anisotropic stresses. 

We assume that the defects create a small perturbation 
on a homogenous and isotropic cosmological background. 
The metric may be written as a small departure from the 
Friedmann-Lemaitre-Robertson- Walker (FLRW) form 
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with dx^ = dt the conformal time and a{t) the scale fac- 
tor. GWs are represented by the transverse and traceless 
(TT) parts of the metric perturbations fijj^ , satisfying 
dihj^ = hj^ — 0. Expanding the Einstein equations to 
first order in hj^ ^ we obtain 
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where hj!^ = ahj^ , and Uj^ is the TT part of Tij. 
The spectrum of energy density of a stochastic GW back- 
ground in comoving momentum k is given by 
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with |/ifc(f)p the power spectrum of hj^ . The solution 
in Fourier space to Eq. ([3|) is 
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with ti the initial time of GW production, i.e. hlj (ti) — 
hj^^{ti) = 0, and Q{k,t,t') the retarded Green's function 
associated to the differential operator in the left hand 
side of Eq. ([3]). At sub-horizon scales {kt,kt' ^ 1), 
Q{k,t,t') = fc^^sin(A:(t — t')). Averaging over a time 
St = 2TT/k, the GW spectrum becomes 
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X coB{k{t' -t"))n^{k,t',t"), (6) 
where H^ is the unequal time correlator (UETC) of H^""", 

(H;^T(k, OHT^tk', t')) = (27r)^H2 (fc, t, t')Sn (k - k') . (7) 

The correlator H^(fc,ii,i2) can be obtained in general 
from field theory simulations. If the source is scaling, 
then the UETC can only depend on k through the vari- 
ables xi — kti and X2 = kt2. From dimensional analysis 
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with the factor 4 chosen so that agrees with the tensor 
UETC of Ref. [l^- Using this form of the correlator 
and the fact that in a radiation background the scale 
factor normalized to unity today can be written as a{t) = 
y/TT^Ho t, we obtain at sub- horizon scales x = kt ^ 1, 
that the spectrum of GW becomes 



dp. 



GW 



{x,t) = rirad- 



4 MlH^ 



F^ix) 



(9) 



dlogfc'"'''^' ' TT a{t)^ \MpJ 
F'^{x) = 2 dxi / dx2y/xiX2 cos(a;i — X2) C'^{xi,X2) (10) 



At subhorizon scales C'^ {xi,X2) is peaked near xi = 
X2, and decays along the diagonal as a power law (see 
e.g. 0)- It also decays away from the diagonal due to the 
lack of time coherence of the source [l8| . Hence the con- 
vergence of the integration is guaranteed for fast enough 
decays. That implies that F'^ {x) becomes more and more 
insensitive to its upper bound of integration, approach- 
ing asymptotically a constant value for x ^ 1. In other 
words, F^ — F^{x 00) is a constant. As a conse- 
quence of this, the GW spectrum at subhorizon scales 
becomes scale-invariant. 

For every type of defect there is indeed a function 
C"^(a;i, 2:2), and thus a well-determined value F^, which 
characterizes the amplitude of the GW background. In 
particular, redshifting the amplitude today and using 
SH^Mp = 8npc, we obtain 



(4) i^Gw(fc) = — 
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That is, the background of GWs produced during the 
radiation era by the evolution of any network of defects in 
scaling regime, is exactly scale-invariant. The amplitude 
is supressed by the fraction Jlrad, and it depends on the 
vev as (v/Mp)'^, and on the shape of the UETC, which 
ultimately modulates the amphtude through F^. We 
can identify the vahie of in the large N analytical 
calculation ^ as (32/ 3) F^{N) ~ 650/N. 

The method of calculating UETCs from lattice field 
theory simulations is well-documented [l3, [l9l - [2l| . We 
consider a model with a global 0{N) symmetry, spon- 
taneously broken to 0(A^ — 1) in the ground state. We 
take a scalar field $ with TV (real) components, $ = 
((/?i, ipn)'^ /V^- The lagrangian of the model and the 
energy-momentum tensor are given by 

£($) A($T$_„2/2)^ (12) 

T^,(x, t) = 2d^<^>^d,<i> - g^, £($) (13) 

where = 1/2^^ i^^. After symmetry breaking the 
scalar is very close to its vacuum expectation value. 

Taking the spatial Fourier transform of T^, the two 
tensor polarizations [A = 1,2) contributing to the GW 
source are defined as 

Sl{Kt)^^J2M,jT,,ik,t), (14) 

where the projectors obey J^a ^ij^hn = ^ij,im = 

PiiPjm - ^f^PijPinii with Pij = 5.1J - kSj- Here Kj,im 
is the projector onto the TT part of . The UETC is 

C^{x,,x,)^\Y,{Sl{^^M)Sl{Kt2r), (15) 

^ A 

where the average is taken over a set of numerical sim- 
ulations and a shell in Fourier space. The UETC obeys 
the symmetry C'^ {xi,X2) — C'^ {x2,xi). 

The algorithm for the numerical simulations solves 
the Klein-Gordon field equations obtained from the la- 
grangian above, on a periodic Cartesian grid and using 
a 7-point stencil for the 3D-Laplacian and a leapfrog for 
the time evolution. We keep a^X constant to maintain 
a constant comoving scalar mass m = VXv, with A = 1 
and w = 1 The grid size was 1024^, the spac- 

ing was Aa; = 0.5, and the timestep At = 0.2Ax. The 
fields are typically initialised at conformal time t = 1 
with independent random values constrained to lie on 
the {N - l)-sphere <P^^ = and with $ = 0. The 

system is evolved initially with a period of variable diffu- 
sion or dissipation, so that it relaxes quickly to scaling. 
The UETCs are constructed by multiplying the Fourier 
transforms from 50 logarithmically spaced times in the 
range t^ct < < ^cnd, with the one taken at a refer- 
ence time t2 = ircf- For N > 2, (trefj^end) = (64,232). 
Strings need more time to reach scaling, so for TV = 2, 
{tref,tend) = (150,300). To obtain the power spectra we 
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FIG. 1: ETCs from simulations, NE'^^^^{x) (solid black), eval- 
uated at t = 64 {t — 150 for N — 2). From highest to lowest, 
we give iV = 2, 3, 4, 8, 12 and 20. Lowest of all is NE'^^{x) 
(dashed). For A/' = 4 we also show the ETCs at f = 232, 
and times in between in grey, to demonstrate the excellent 
scaling (for N = 2, due to later onset of scaling, we plot from 
X > 1.84). The error bars on the numerical curves give the 
la variation over all runs, and are barely visible. 



average over a shell of width Ak — 2tt/L, where L is the 
side length of the simulation volume. 

To compare the dynamics in our lattice with the large 
N analytical calculation, we computed the equal time 
correlator (ETC) E'^{x) = C'^{x,x) for each of TV = 2, 
3, 4, 8, 12 and 20. The theoretical ETC in the large TV 
limit scales as £^{^(a;) oc 1/TV. Thus, in Fig. [T] we plot 
NE^^^{x) at irof, obtained (averaging over 20 realiza- 
tions) for each TV, with NEJ^{x) for comparison. 

Computing the numerical to theoretical ratio of ETCs 
at the scale x = it when one full wavelength enters the 
horizon, Tjv = E'^^^-^{'k) / EJ^{'k) , we find the values ex- 
hibited in Table HI The numerical ETC approaches the 
theoretical prediction as TV grows, as shown by the ap- 
proach of Tat to unity as TV increases. For the case of 
cosmic strings (TV = 2) the numerical ETC is a factor 
^ 40 bigger than the theoretical one, signalling the break- 
down of the large TV approximation. We do not expect 
the large TV approximation to apply for N — 2: not only 
is TV small, but there are string defects which invalidate 
the mean-field analysis from the start. The reason that 
the scaling density is so much larger than the large TV 
value is however unclear to us. 

In Fig. [2] we compare the GW amplitude today JIqw 
given by the large TV theoretical calculation Eq. ([T]), ver- 
sus the amplitude f^^™ obtained from Eq. pd|) , with F^ 
calculated from the numerical UETC for each scenario. 
We see that rJo^/r^Q^ approaches unity as TV grows, 
showing that the numerical amplitudes converge to the 
large TV analytical result. Today's amplitudes for each 
different TV considered are summarized in Table ID If we 
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FIG. 2: Ratio of the numerical GW amplitude rJow' to the 
large analytical calculation fiav^ (see Table and a fit to 
1.1+45/iV^. The error bars give the la variation over all runs. 

parametrize the numerical amplitude for > 4 as 

f^GW = ^'^w («o + ^ + ^ + •■•) , (16) 

we find a good approximation with oq — 1.1, 02 — 45 
(see Fig. [2]) with negligible oi. We believe that the 10% 
deviation from unity of oq is due to a systematic uncer- 
tainty in our numerics, most likely a finite volume effect. 
Thus, the numerical GW amplitude approaches the large 
N result faster than naively expected, as cx albeit 
with a large coefficient. The convergence refiects the be- 
haviour of the overall scale of the UETCs, as measured by 
Ttv, although we see some TV-dependence in the UETC 
width, which we shall report on in the future. Strings are 
well above this trend, by a factor of about 100. We can- 
not be more precise at this stage, as there is a systematic 
uncertainty arising from the extrapolation of the UETC 
to a; = 0. We estimate this to be of order 50%. 



iV 


2 3 4 8 12 20 


Tiv 


36 4.5 3.1 1.7 1.4 1.3 


(-)num /(-)th 
"GW/^'GW 


130 7.3 3.9 1.8 1.4 1.3 



TABLE I: Values of the numerical ETCs at x = vr, and GW 
amplitudes today, normalized to the large A'' calculation. The 
fluctuation in the amplitudes over the 20 realizations is less 
than 10%, except for Af = 2 where it is ~ 20%. 

In this letter we have clarified the origin of the scale- 
invariance of the GW background calculated in [l3, [ij] 



for the case of non-topological global textures. More im- 
portantly, we have generalized the result: a scale invari- 
ant background of GW is expected from any scaling cos- 
mological source during the radiation era. In particular, 
global defects, independent of their topology, and cosmic 
gauged strings (local or semi- local), enter into a scaling 
regime, and produce a scale-invariant (i.e. frequency- 
independent) GW power spectrum according to Eq. , 
whose amplitude depends on the defect type. 

We performed numerical simulations of the self- 
ordering dynamics of an 0{N) scalar field, showing that 
the GW power spectrum approaches the large N predic- 
tion at a rate consistent with iV^^ (with a surprisingly big 
coefficient). For example, for = 4 the GW power spec- 
trum is approximately four times larger than the large A'^ 
prediction. For strings, the factor is of order 100. 

We note that global strings {N — 2) decay by emis- 
sion of massless [2^ and massive scalar radiation, both 
from infinite strings and loops, at a rate proportional 
to [v/Mpf. Hence the GW emission, whose power is 
proportional to {v/Mp)'^, is not a significant source of 
energy loss. Global strings therefore do not behave like 
local strings in the Nambu-Goto approximation, which 
decay into GWs alone, via emission from sub-horizon 
size string loops. The amplitude of this background de- 
pends sensitively on the as-yet uncertain loop size dis- 
tribution (see [2j, [i^l for a recent well-referenced inves- 
tigation). We emphasise that the background we pre- 
dict arises from long strings and short-lived horizon-size 
loops, and has not been considered before. While sub- 
dominant for Nambu-Goto strings, it forms an irreducible 
minimum for strings decaying by particle emission. 

It will be interesting to calculate the GW power spec- 
trum from gauge cosmic strings, where numerical simula- 
tions show that the ETC decays much more slowly. The 
GWs can contribute appreciably to the relativistic en- 
ergy density, with important implications for the cosmic 
microwave background power spectrum [25| . 
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